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1.0  INTRODUCTION 

The  computation  of  the  flow  field  with  viscous  effects  over  an  airfoil  or  similar  body 
is  usually  carried  out  by  an  iterative  method.  The  first  step  is  to  obtain  an  inviscid  outer 
solution  of  the  basic  profile.  A  displacement  thickness  is  then  calculated  by  using  the 
resulting  velocity  or  pressure  distribution  in  a  boundary-layer  solution.  This  displacement 
thickness  is  added  to  the  initial  body  profile  to  form  an  equivalent  body  and  the  process 
is  repeated.  Convergence  is  normally  achieved  with  a  few  iterations. 

This  iterative  procedure  is  not  often  successful,  however,  if  significant  flow  separation 
occurs.  Standard  boundary-layer  methods,  if  used  in  a  separated  flow  region,  produce  highly 
inaccurate  displacement  thicknesses.  This  is  particularly  true  in  transonic  flows  in  which 
a  shock  appears,  due  to  the  strong  pressure  gradients  induced  by  the  shock  and  possible 
shock-induced  separation. 

Two  examples  of  transonic  flows  exhibiting  the  shock  boundary-layer  interaction  are 
shown  in  Fig.  1 .  The  first  case,  a  modified  circular  arc  bump  mounted  on  a  wind  tunnel 
floor,  was  investigated  by  Alber.  et  al.  (Ref.  1 ).  The  second  velocity  distribution  consists 
of  data  from  the  test  of  an  airfoil-shaped  bump  mounted  on  the  floor  of  AEDC 
Aerodynamic  Wind  Tunnel  (IT),  a  test  performed  by  F.  L.  Heltsley  in  May  1974.  While 


X,  feel 


a.  Circular  arc  bump 

Figure  1.  Comparison  of  experimental  velocity  distributions 
with  velocities  obtained  by  inviscid  computation. 
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b.  Airfoil-shaped  bump 
Figure  1.  Concluded. 

this  second  flow  is  not  fully  separated,  it  presents  the  same  difficulties  due  to  the  greatly 
increased  boundary-layer  thickness  downstream  of  the  shock. 

A  comparison  of  the  experimentally  measured  velocities  with  those  obtained  using 
Murman's  inviscid  solution  (Ref.  2)  clearly  indicates  the  error  resulting  from  ignoring  the 
shock  boundary-layer  interaction.  The  large  increase  in  boundary-layer  thickness  caused 
by  the  shock  forms  an  equivalent  body  greatly  different  from  the  original  body. 

As  part  of  an  effort  to  extend  an  iterative  method  of  solution  to  separated  flows 
of  the  type  illustrated  in  Fig.  1,  a  computer  program  for  a  turbulent,  compressible, 
boundary-layer  method,  capable  of  carrying  out  computations  in  a  separated  flow  region, 
has  been  developed  and  tested.  The  procedure  used  for  extension  through  the  separation 
point  follows  that  used  by  Kuhn  and  Nielsen  (Ref.  3).  Instead  of  specifying  a  pressure 
or  external  velocity  distribution,  friction  velocity  is  specified,  so  that  the  external  velocity 
is  obtained  as  a  result.  The  prescribed  friction  velocity  is  adjusted  so  that  the  external 
velocity  computed  by  the  boundary-layer  method  matches  a  velocity  measured 
experimentally  or  an  inviscid  computation.  Alternately,  in  the  present  method,  the 
boundary-layer  thickness  can  be  specified  in  the  same  manner. 

By  prescribing  either  the  friction  velocity  or  the  boundary-layer  thickness,  the 
singularity  occurring  as  skin  friction  goes  to  zero  is  avoided,  more  accurate  boundary-layer 
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velocity  profiles  are  produced,  and  the  displacement  thicknesses  obtained  compare  well 
with  experimental  results. 


2.0  BOUNDARY-LAYER  METHOD  FOR  ATTACHED  FLOW 

The  turbulent,  compressible,  boundary-layer  method  used  is  a  variation  of  the  method 
of  Nash  and  Hicks  (Ref.  4)  using  a  modification  of  Cole's  wall-wake  velocity  representation 
(Ref.  5)  to  replace  the  shear  stress  equation  proposed  by  Nash.  The  Stewartson 
Transformations  (Ref.  6)  are  included  for  compressibility  effects.  This  method  was  modified 
by  Kuhn  for  boundary-layer  analysis  in  a  region  of  separated  flow  (Ref.  3). 

The  reasons  for  choosing  this  method  are  speed,  flexibility,  and  the  previous 
experience  gained  by  Kuhn  in  and  near  separated  regions.  The  speed  is  particularly 
important  as  the  future  use  of  this  method  is  part  of  an  iterative  solution. 

2.1  DEVELOPMENT  OF  THE  BOUNDARY-LAYER  EQUATIONS 

Starting  with  the  continuity,  momentum  and  energy  equations  for  a  turbulent, 
compressible  boundary  layer  (Ref.  3) 

9  9 

(pu)  +  —  (pv)  =  0  (1) 


(2) 


(3) 


where 


S  =  T/T,  -  1 


(4) 


The  Stewartson  Transformations  (Ref.  6)  are  applied  in  the  following  form 
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With  the  assumptions  that  the  fluid  is  a  perfect  gas,  the  laminar  and  turbulent  Prandtl 
numbers  are  unity,  viscosity  is  linear  with  temperature,  and  the  wall  is  adiabatic,  the 
set  of  Eqs.  (1  through  41  can  be  transformed  to 

Ux  ±  Vy  =  0  (9) 


UUx  +  VUy  -  (S  +  1)  UeUc^  -  »^(i3Uy)y  =  0  (10) 


S  = 


Taw 

-= - 1  =  constant 


(11) 


Equations  (9  through  11)  are  combined  using  the  usual  boundary-layer  integral  approach 
(Ref.  4). 

/|l^Ux  -  Uy  l'  UxdT,  -  UeUe,  -  m  Uy)y  |  y" dy  =  0  (12) 

Equation  (12)  is  the  integral  of  the  momentum  across  the  boundary  layer  for  n  =  0,  and 
the  moment  of  momentum  for  n  =  1. 

With  the  addition  of  a  boundary-layer  velocity  profile  representation  and  an  expression 
for  eddy  viscosity,  a  closed  set  of  equations  can  be  formed.  Crocco's  theorem  with  a 
relaxation  factor  of  0.89  is  used  to  determine  density  gradient  across  the  boundary  layer. 
The  velocity  profile  is  taken  in  the  form  given  by  Coles  (Ref.  5)  with  an  exponential 
term  added  for  the  viscous  sublayer: 

U  =  ut[2.5  fin  (1  +  y+)  +  5.1  -  (3.39y+  +  5.1)e-0-37y-"]  +  ^  .  cos  (13) 

where 


y+  =  \\ir\ylv 


(14) 


and 
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Ut 


=  /l^w 
'  Tw  I  \  p  y 


(15) 


The  friction  velocity,  uj,  is  defined  in  this  way  to  allow  for  reverse  flow  in  a  separated 
region.  The  wake  velocity.  UjS,  can  be  eliminated  by  setting  y  =  3  in  Eq.  (13): 


Uj3  =  Uc  -  Ut(2.5  fin  (1  +  5+)  +  5.1) 


(16) 


For  the  unseparated  flow,  the  eddy  viscosities  are  expressed  in  two  forms,  for  an 
inner  and  outer  region  (Refs.  7  and  8). 


ft  =  1  +  0.0533  jcO-41  -  j^l  +  0.41U+  +  y  (0.41U+)2j| 


where 


and 


U+  =  U/ur 


/3o  =  1  +  K  ^1  +  5.5  (-f-y  )*  ^cS*lv 
with  K  taken  as  0.0168  for  favorable  pressure  gradients  and 


(17) 


(18) 


(19) 


K  =  0.013  +  0.0038  exp  ^-6*  -^/15t,v^ 


(20) 


for  dp/dx  positive  (Ref.  3). 

In  a  separated  flow  region,  following  Alber  (Ref.  9),  the  eddy  viscosity  is  given  by 


0  -  0.013  [l  .  5.5  ‘  -  .iL)dy 

yu=o 


(21) 


Substituting  the  velocity,  Eq.  (13),  and  its  derivatives  into  Eq.  (12)  results  in  the 
two  equations 


A]  1  uf  +  A]2  6'  +  A]  3  Uc  =  -Ut|ut|  (22) 

A21  Ur  +  A22  5'  +  A23  E'c  =  -V  /*/3Uydy  (23) 

The  prime  denotes  differentiation  with  respect  to  x.  The  Ajj  coefficients  are  functions 
of  Uj,  Uc  and  6  as  defined  in  Appendix  A. 
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2.2  METHOD  OF  SOLUTION 

Given  the  inviscid  velocity  distribution,  Eqs.  (22)  and  (23)  form  an  initial  value 
problem  in  the  dependent  variables  Ur  and  6.  In  addition  to  the  velocity,  initial  values 
of  U7.  6.  and  the  free-stream  thermodynamic  state  are  required.  Input  is  in  terms  of  the 
physical  variables  and  both  compressible  and  incompressible  results  are  output. 

Equations  (22)  and  (23)  were  integrated  using  a  fourth  order  Runge-Kutta  method 
(Ref.  10)  with  the  step  size  in  x  of  the  order  of  5.  The  interval  for  integration  across 
the  boundary  layer  was  6/20  for  the  cases  shown.  Accuracy  is  not  very  sensitive  to  the 
step  size  chosen  if  it  is  in  that  range.  The  time  required  for  each  streamwise  integration 
step  is  approximately  0.05  sec  on  an  IBM  370/165.  A  listing  of  the  computer  program 
is  given  as  Appendix  B. 

Note  that  the  problem  is  formally  the  same  if  either  U7  or  6  is  specified  and  the 
external  velocity  is  treated  as  a  dependent  variable. 

2.3  TESTING  OF  THE  VALIDITY  OF  THE  METHOD 

The  performance  of  the  present  boundary-layer  method  was  checked  by  comparison 
with  other  boundary-layer  methods  and  with  experimental  data.  Figure  2a  shows  the 
pressure  distribution  on  the  upper  surface  of  a  6-in.  C-141  airfoil  taken  in  Tunnel  16T 
(Ref.  1 1 ).  The  boundary-layer  characteristics  shown  in  Figs.  2b  through  e  were  calculated, 
using  the  C-141  data,  by  the  present  method,  the  method  of  Nash  and  Hicks  (Ref.  4), 
the  method  of  Patankar  and  Spalding  (Ref.  12)  as  modified  by  High  and  Felderman  (Ref. 
13),  and  by  Adams  (Ref.  14). 

There  is  considerable  difference  in  values  for  some  of  the  results,  particularly  shape 
factor  and  friction  coefficient  (Figs.  2b  and  c).  This  difference,  with  the  exception  of 
shape  factor  variation,  is  not  unusually  large,  even  for  comparison  of  incompressible 
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methods  in  moderate  pressure  gradients  (Ref.  15).  The  momentum  thicknesses  compared 
in  Fig.  2d  are  in  better  agreement,  and  the  displacement  thicknesses  (Fig.  2e)  are  quite 
close.  Since  the  result  of  principal  interest  for  present  purposes  is  displacement  thickness, 
the  agreement  is  considered  satisfactory. 

The  comparison  with  experimental  data  is  shown  in  Figs.  3a  and  b.  The  computations 
are  based  on  the  experimental  velocity  distribution  over  the  circular  arc  shown  in  Fig. 
la.  The  experimental  data  were  presented  by  Alber,  et  al.  (Ref.  1).  The  agreement  with 
experiment  for  both  the  friction  velocity  and  the  displacement  thickness  is  quite  good 
up  to  the  strong  pressure  gradient  near  the  shock.  Then  the  accuracy  becomes  poor, 
especially  after  separation. 


CHORD  STATION,  X/C 

a.  Pressure  coefficient  on  the  airfoil 

Figure  2.  Comparison  of  results  from  four  boundary-layer  methods  using 
the  pressure  distribution  on  a  C-141  airfoil. 
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X,  feet 

b.  Displacement  Thickness 

Figure  3.  Boundary-layer  characteristics  computed  from  the  experimental 
velocity  distribution  on  the  circular  arc  bump. 

3.0  PROCEDURE  IN  A  SEPARATED  FLOW  REGION 

Three  problems  are  most  evident  in  boundary-layer  computations  involving  a  shock 
and  separation.  The  first  is  the  strong  pressure  gradients  in  the  vicinity  of  the  shock. 
The  boundary-layer  equations  are  not  developed  for  use  in  strong  pressure  gradients  and 
the  assumptions  made  in  the  derivations  are  not  very  appropriate.  The  second  problem 
is  prediction  of  the  point  of  separation.  Boundary-layer  methods  typically  do  not  predict 
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separation  accurately,  as  the  equations  are  singular  at  that  point  (Ref.  16),  and  while 
some  prediction  schemes  such  as  Alber’s  (Ref.  9)  and  Stratford's  (Ref.  17)  are  better, 
the  accuracy  in  the  transonic  range  is  still  inadequate.  Finally,  computing  boundary-layer 
characteristics  in  a  separated  region  with  velocity  specified  normally  results  in  very  poor 
accuracy.  This  third  problem  is  treated  here  by  specifying  either  the  friction  velocity  or 
the  boundary-layer  thickness  in  Eqs.  (22)  and  (23)  and  solving  for  external  velocity  as 
a  dependent  variable. 

3.1  METHOD  OF  COMPUTATION 

Since  Eqs.  (22)  and  (23)  are  in  terms  of  three  unknowns,  external  velocity, 
boundary-layer  thickness,  and  friction  velocity,  one  of  these  variables  must  be  prescribed. 
Usually  the  external  velocity,  either  measured  experimentally  or  from  an  inviscid  solution, 
is  specified  and  6  and  Ut  are  calculated.  Near  separation,  this  results  in  large  inaccuracies. 
Typically,  the  skin-friction  coefficient  (or  friction  velocity)  decreases  in  value  but  does 
not  reach  zero  (Ref.  16).  Similarly,  other  boundary-layer  parameters  fail  to  reach  accurate 
values.  Figure  3  is  an  example  of  this  behavior  as  exhibited  by  the  present  method. 

The  accuracy  of  the  results  can  be  improved  by  specifying  friction  velocity  or 
boundary-layer  thickness  and  solving  for  external  velocity  as  a  dependent  variable.  In  order 
to  produce  the  desired  experimental  or  calculated  (inviscid)  external  velocity,  a  trial  and 
error  procedure  is  necessary.  The  specified  friction  velocity  or  boundary-layer  thickness 
is  varied  until  the  resulting  external  velocity  matches  that  desired. 

Accuracy  is  improved  both  by  forcing  the  boundary-layer  velocity  profiles  to  assume 
a  more  correct  shape  and  by  avoiding  the  singular  behavior  which  occurs  as  skin  friction 
approaches  zero  when  external  velocity  is  specified. 

3.2  RESULTS 

Figure  4  illustrates  a  case  in  which  the  friction  velocity  is  specified.  The  body  is 
the  circular  arc  profile  shown  in  Fig.  la.  Figure  4a  gives  the  specified  u,-  as  compared 
with  experimental  values.  This  Ut  distribution  is  close  to  the  experimentally  determined 
value,  and  the  computed  external  velocity  is  also  near  the  experimental  value  (Fig.  4b). 
From  Fig.  4c  it  can  be  seen  that  the  displacement  thickness  is  determined  with  good 
accuracy. 

Figure  5  illustrates  the  same  process  with  boundary-layer  thickness  specified.  Figure 
5a  is  the  prescribed  boundary-layer  thickness.  Fig.  5b  is  the  resulting  external  velocity, 
and  Fig.  5c  compares  the  computed  displacement  thickness  with  experimental  values.  Again, 
the  agreement  is  quite  close. 
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o  ALBER,ET  AL.,REF.  I 
-  PRESENT  RESULTS 


a.  Specified  friction  velocity 


b.  External  velocity 


c.  Displacement  thickness 

Figure  4.  Boundary-layer  characteristics  computed  for  the  circular 
arc  bump  with  friction  velocity  specified. 
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o  ALBER,  ETAL.,REF.  I 
-  PRESENT  RESULTS 


a.  Specified  boundary-layer  thickness 


b.  External  velocity 


X,  feet 

c.  Displacement  thickness 

Figure  5.  Boundary-layer  characteristics  computed  for  the  circular 
arc  bump  with  boundary-layer  thickness  specified. 

Figure  6  compares  computed  and  experimental  boundary-layer  velocity  profiles  for 
the  flow  when  boundary-layer  thickness  is  specified.  Figure  6a  is  a  profile  upstream  of 
separation,  Fig.  6b  is  in  the  separated  How,  and  Fig.  6c  is  after  reattachment.  Similar 
accuracy  is  obtained  when  friction  velocity  is  specified. 
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—  PRESENT  RESULTS 


Figure  6.  Computed  and  experimental  boundary -layer  velocity  profiles 
for  the  circular  arc  bump. 
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Figure  7.  however,  shows  a  difficulty  that  can  arise  while  using  this  method.  The 
computations  for  this  Figure  are  based  on  the  experimental  velocity  distribution  shown 
in  Fig.  lb.  Figures  7a  and  b  demonstrate  a  specified  friction  velocity  and  the  computed 
external  velocity.  Similarly.  Figs.  7c  and  d  give  a  specified  boundary-layer  thickness  and 
the  associated  velocity.  Both  calculated  velocities  agree  with  experimental  data;  however, 
Fig.  7c  shows  the  discrepancy  between  the  displacement  thickness  resulting  from  these 
two  computations,  hven  though  the  practical  difference  between  these  displacement 
thicknesses  is  not  large.  Fig.  7e  indicates  an  ambiguity  that  can  occur  in  these  solutions. 


O  ALBER,  ETAL.,REF.I 
-  PRESENT  RESULTS 


a.  Specified  friction  velocity 


b.  Resulting  velocity  profile 

Figure  7.  Calculation  of  boundary-layer  characteristics  for  a  flow  near 
separation  over  an  airfoil-shaped  bump. 
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O  ALBER,ETAL..REF.I 

-  PRESENT  RESULTS 
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c.  Specified  boundary-layer  thickness 


d.  Resulting  velocity  profile 
Figure  7.  Continued. 
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e.  Comparison  of  displacement  thicknesses  computed  by  this  method 

Figure  7.  Concluded. 

4.0  CONCLUSIONS 

The  purpose  of  this  work  was  to  develop  a  computer  code  for  a  compressible, 
turbulent,  boundary-layer  method  which  could  be  used  for  computations  in  a  region  of 
separated  flow. 

The  boundary-layer  method  presented  here  is  fast  and  relatively  accurate  for  attached 
flows  in  moderate  pressure  gradients.  As  with  most  available  techniques,  the  accuracy  is 
poor  in  very  strong  pressure  gradients,  such  as  those  in  the  vicinity  of  a  shock. 

Following  Kuhn  and  Nielsen,  the  solutions  can  be  extended  to  separated  flow  regions 
by  specification  of  friction  velocity  as  the  independent  variable  in  the  boundary-layer 
equations  and  solving  for  external  velocity  as  a  dependent  variable.  In  the  present  program, 
the  boundary-layer  thickness  can  also  be  specified  instead  of  friction  velocity.  When  either 
friction  velocity  or  boundary-layer  thickness  is  specified  so  that  the  resulting  external 
velocity  matches  an  experimental  or  calculated  inviscid  external  velocity,  accurate 
displacement  thicknesses  can  be  obtained. 
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APPENDIX  A 

DERIVATION  OF  DIFFERENTIAL  EQUATIONS 

Using  Eqs.  (12)  and  (13)  one  obtains: 

f  juUx  -  Uy  ^  Ux  dTj  -  (S  +  1)  Ue  Ue^  -  V  (/3Uy)yf  yndy  = 
0  0 

0  (A-1) 

n  =  0,  1 

and 

U  =  Ur  [2.5  8n  (1  +  y+)  +  5.1  -  (3.39  y+  +  5.1)  e-0-37y-"]  +  ^1 

7ry\ 

-  cos  —j 

where 

(A-2) 

|rw|  (  p  ) 

(A-3) 

Up  =  Ue  -  Ut  (2.5  fin  (1  +  6+)  +  5.1) 

(A^) 

3 

II 

(A-5) 

8+  J“'l* 

V 

{A-6) 

Then  let 

Fl(y)  =  2.5  fin  (1  +  y+)  +  5.1  -  (3.39  y+  +  5.1)  e-0-37y* 

(A-7) 

2.5 

F2(y)  =  - - —  +  (1.254  y+  -  1.502)  e-0-37y+ 

1  +  y"^ 

(A-8) 

Y  (>  ■  (v)) 

(A-9) 

F4(y)  =  -^  sin  -g- 

(A-10) 
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SO  that  Eq.  (A-2)  becomes 

U  =  urFiCy)  +  UjjFaCy)  (A-11) 

Taking  the  derivatives  of  Eq.  (A-11)  with  the  additional  definition 

FsCy)  =  Fi(y)  +  y+  F2(y)  (A-12) 


Ux  = 


uf  (FsCy)  -  FsCy)  F5(6))  -  S'  F4(y)  +  Ut  FsCy)  F5(y)^+  Ui  FjCy) 

(A-13) 

Uy  =  Ut  F2(y)  + 


5 


F4(y) 


(A-14) 


Defining 

Fs  =  Fsfy)  -  F5(6)F3{y)  (A-15) 

F7  =  -  F4(y)  -  u^F2(6)F3(y)  (A-16) 

Equation  (A-1)  can  then  be  written 

fi  /  y 

/  |u[u^F6  +u;  F3(y)+  S’F^]  +  Uy  J  [u^  +  U'  F3(y)  +  6'F7]  d»? 

0  ^  0 

-  (s  +  DUeUe  -  V  (/3  Uy)y|y"dy  =  0  (A-17) 

Then  the  coefficients  for  Eqs.  (19)  and  (20)  become 


Aki  =  J 

/uF6  +  Uy/Fs  d77^ 

y(k-l)  dy 

(A-18) 

0 

V  0  / 

r® 

Ak2  =  j 

(uF7  +  Uy  f ^  F7  dr?) 

y(k-l )  dy 

(A-1 9) 

0 

V  •'o  / 

.6 

/  -y 

\ 

Ak3  =J 

UF3(y)  +  Uy  1  F3(y)dT7  -  (s  +  1)  UeUMyt'^  n  dy 

(A-20) 

0 

N  0 

/ 

These  expressions  were  evaluated  numerically. 
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APPENDIX  B 
FORTRAN  IV  PROGRAM 

B-1  DESCRIPTION  OF  PROGRAM  FUNCTIONS 

Given  initial  values  and  the  external  velocity,  this  program  computes  characteristics 
of  a  turbulent,  compressible  boundary  layer.  For  a  given  region  in  the  flow,  either  friction 
velocity  or  boundary-layer  thickness  can  be  specified  and  external  velocity  computed  as 
a  dependent  variable. 

B-2  INPUT 


All  input  except  specification  of  friction  velocity  or  boundary-layer  thickness  is  read 
-  in  Subroutine  Input.  The  specified  variables  are  read  in  Subroutine  Loop. 

Card  Variables  Format  Description 

1  M,  N  (315)  M  is  the  number  of  streamwise  grid 

points,  corresponding  to  the 
number  of  external  velocities 
input.  N  is  the  number  of  points 
'  taken  across  the  boundary  layer. 

N  =  21  is  recommended. 


2 


3 


LS,LD,LT 


C,PI,RHOI, 

UI,DLC(1) 

UM,XO, 

UTAU(l) 


(315)  LS  is  the  first  streamwise  grid 

point  at  which  friction  velocity 
is  specified.  LD  is  the  last 
point  specified.  If  LS  =  LD  then 
boundary-layer  thickness  is 
specified  and  LT  is  the  last 
point.  If  LS  is  greater  than  M 
neither  variable  is  specified. 

(8F10.0)  C  is  the  length  of  the  computed 
area  in  feet.  C/(M-1)  is  the 
streamwise  computational  interval. 
PI,  RHOI  and  UI  are  the  reference 
free-stream  pressure,  density 
and  velocity  in  Ib/sq  ft,  slugs/cu 
ft  and  ft/sec.  DLC(l)  is  the 
initial  boundary-layer  thickness 
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Card  Variables  Format 

3  (Continued) 


4-K  UEC(I)  (7F10.0) 

(K+D-  UTAU(I)  (7F10.0) 

END  or 

DEL(I) 


Description 

in  feet.  UM  is  the  reference 
absolute  viscosity  in  lb  sec/sq 
ft.  XO  is  the  stream  wise  location 
of  the  starting  point.  UTAU(l) 
is  the  initial  value  of  the 
friction  velocity,  V^w/p,  in 
ft/sec. 

External  velocity  in  ft/sec.  Must 
be  given  for  each  X  station,  but 
dummy  values  may  be  used  if  UTAU 
or  DEL  is  specified  at  that 
location. 

Either  friction  velocity  (ft/sec) 
or  incompressible  boundary-layer 
thickness  (ft)  as  determined  by 
card  2.  Input  values  for  alternate 
points;  the  program  will 
interpolate. 


B-3  OUTPUT 

Output  is  in  the  same  units  as  input.  Print  interval  is  controlled  for  variables  K4 
and  K6  in  Input.  On  the  first  page  are  starting  conditions  and  initial  profiles,  in  compressible 
and  incompressible  form. 

The  next  block  gives  results  at  the  desipiated  streamwise  locations.  The  following 
variables  are  printed; 

N  Number  of  the  station 

XC  Compressible  (physical)  streamwise  location 

UC  Compressible  external  velodty  input,  or  computed  if  UTAU  or  DEL  is 

specified 

DU/DXC  Derivative  of  the  input  external  velocity,  does  not  correspond  to  UC  in 
regions  where  UTAU  or  DEL  is  specified 

DELC  Compressible  boundary-layer  thickness 
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DEL* 

Compressible  displacement  thickness 

TH 

Compressible  momentum  thickness 

H 

Shape  factor 

CF 

Skin-friction  coefficient 

NU 

Kinematic  viscosity 

X 

Transformed  incompressible  streamwise  location 

UEI 

Incompressible  external  velocity,  computed  value  if  UTAU 
specified 

or  DELI  is 

DU/DXI 

Incompressible  velocity  derivative,  also  computed  if  UTAU 
specified 

1 

or  DELI  is 

DELI 

Incompressible  boundary-layer  thickness,  this  is  one  of  the  variables  that 
can  be  specified 

UTAU 

Friction  velocity  -  this  is  the  other  variable 

specified.  1  I  ^ 

that  can  be 

UBET 

Wake  velocity  from  Eq.  (16) 

TW 

Wall  temperature 

The  last  print  block  gives  compressible  velocity  profiles  in  the  boundary  layer.  The 
profiles  are  given  at  streamwise  locations  corresponding  to  those  having  the  same  number 
in  the  preceding  print  block.  Values  printed  are  U/Ue  in  thousandths,  with  the  decimal 
and  leading  zeros  omitted.  The  spacing  across  the  boundary  layer  is  constant  in  compressible 
coordinates. 
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i_ 

TMlb  PriObHAr*  CO«PultS  A  I  lIKrtiJLLli  T  nOilnOAHY  LAYtH  fir  Tnt  METriOl) 

c 

c 

c 

bIVtN  aY  RIJMN  AnU  imIlLSi-N,  IySILAl)  UP  sHtClPYlNb  Tut  VPLOCITY 

PMOUlM.  IjTAUIPmICI  KIN  YrLuClfY)  UH  UEL I  A  (MUUNOAR  Y  LAYER  THICKNESS) 

CAN  HP.  SHtCIUtU  In  (jPSIbNAltll  RLvjIONS.  INRUT  IS  IN  PttT.  SLUbS 

AND  seconds  0<s  a  Slo^lLAR  CONPAriaLt  SYSTEM. 

c 

c 

input  WAkIAmLES:  m  is  TmL  NUnoLK  up  X  YLLuCITItS  INPUT.  ThtSE 

must  ml  spaced  At  cVPiJ  IuTPMYAlS.  .1  IS  TmL  NOAptn  OE  bKlL)  POINTS 

c 

c 

IN  THE  Y  DlMtCTlON.  POR'-AT  I2JS).  LNltP  LS.LT.LU  (3IS).  LS  IS 

ImE  X  bHIU  POINT  Ml  WHICH  UTAU  IS  P IPST  SPLCIPItU.  LT  IS  ThE  LAST 

c 

JL 

bPIU  POINI  POP  SPLLIPYIhh  UIAU.  ip  LSaLI  THEN  ULLTA  IS  SPtCIEIEI) 

AND  LO  IS  TmE  LaSI  NRIO  POInT,  nlXT  CAPO:  ChOPU  LEnOIm.  PREE 

c 

c 

stream  PRcSSUWt .OEnSI TY  aN|1  VElUCIIY*  initial  bOUNuARY  LAYER 

Thickness.  AbSULuft  viSv.osliY.  Initial  x  station  and  initial  utau 

c 

c 

(PRICTION  VELOCITY).  EOKmaT (bP 1 0 . tf ) .  InPUT  VELOCIIY  AT  EACH  X 
STATION.  FORMA  1  ( /r 10.0) .  EnTPR  (SlARTlNb  NEXT  CARU)  UTAU  OR  OELTa 

c 

c 

IP  APPlICaRI.E.  enter  every  other  point— PRObRAM  xILL  INTERPOLATE. 

FORMAT  1  /P  1 0.0)  . 

c 

MAIM  PRObkAll 

COMMON  YUinl).V(loniljL(10l).UllQil  .KriO(  1  Ull  .  X  T  1 V  0  H  .  X  (  ?n  1  1  . 

1  UEC  (?un  (UE  (<:01 )  .UXC  (col  >  tUX  (x>01 )  .ULC  t^Ul )  «ObL  (^un  •UTAUCi'OU  • 

UTP  tg(il)  iXiM  iPl  .Al»WtiUl  iHbtALiMhULiNtMiCtOXtUYiTMi _ 

c 

3  UELU.Tk.  UhE  f  f  (iNU  .LIn)  .ijn.RR.I  T.I  I.TE.Pi.TS 
establish  CONSIaNIs _ 

KO(ZCt^t) »!♦. 1 /b* (Zf /ML >**?*»!- (/C/Zfe ) **2) 

L (l)sl 

E(£)*2 

jc=o 

1 

KlsO 

V  in  so 

U(1)*0 

Y  r  1 1 1 su. 

UC ( 1 ) BU . 
j=i 

B 

CALL  Input  and  ESTahcIsh  initial  values 

Call  Inpui 

'iAMal.4 

QAM1*GAk-1 

Tl»HR*Pl/RMnI 

AlaSURT  (  1  .4*P1/Pri0l  ) 

Call  ext(i) 

PHOSa  lKH01**bAM  1  PfeAM  1  *  (  HHO  1  **bAM/Pl  )  •  1  LJ I  «»?  1  /  T  ?.  »GAM  )  l»4(l/bAMl  ) _ 

AS  =  St>«r  ♦AlBBg) 

PSsRrtUS*AS**2/bAM. 

TS*WH*MS/MHOb 

AwaSUMT  (r,AM*Ft/kHO(  1)  1 

lJE(l)«UtCin*«X/ME 

UMETaiiEtJ)  -OTmUIJ)  »H  t  DLL  (  J  )  »0T  mU  (  J) /Onu) _ 

L 

ITERATE  POM  initial  VELOCITY  PRuHLE 

UU  11  ia7.N 

Yr(n*YT  <I-1)  ♦!)> 

V(I)aYtl-l)+UY 

U<i  (a'JIttU  1  (mHS(UT  AU(J)  )*Y  ( l)/ONU)  ♦.5«U«tT*(l-CIJS(3.l4l59» 

1  Y(l)/OtL(J)n 

lJC(I)aUlI)*4t/Ai 

1)  MrlOinaKnUt/MOdiLin  tUtCI  n  1 

iSC=N/-i 


14  intSaJClNC) 
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UO  ls^*N 

- 12  YH)»Y(l-l>»(KMOtl)»RnO(l-H  1  *At«>U¥/ IHnUl»Al»?l _ 

UtL(l)*Y(N) 

_ Ubt.T«Ufc.m-»TAUlil«Fl  (UfcL(l)>UTAU(ll/iiNU) _ ^ _ 

OU  13  1>?(N 

- U(11»UTAH(J1«H  (AHSCUTAUtJl  ]*Y  (I)/&NU)».^*LJBtLT»tl-CUSt3.1»mg* 

1  Y(I)/UtL(J)n 

- UCm»U(l)>At./Al _ _ 

13  HhO (I > cKHoe/KU (OL ( 1 > tutc ( 1 ) ) 

_ JC»JC«1 _ _ 

IF  (JC.bT.?n)  60  TU  lb 

- IF  (AMblDTFS-UCtuCn  .6T..m)  L.U  TO  14 _ 

15  OliL(l>«Y(iM) 

- uaE.T«ubm-uTAUti)>Fi  (0tL(ii«LjrAU(n/6NU> _ 

MRlT£(bt200n  JCtOritTtO  [  AIM  1 ) 

_ 00  6  ^ 

Mn0(l)»HHOg/H0(0C«n  »ULC(1»  » 

_ Ytl)»Yll-H*(HnUlH*Hrt01I-n  l*At»OY/(HhOl»Ai*?l _ 

6  U  1 1  )«UT  AU(  J»  *H  I  aHS  (OlAU  ( J)  1  *Y  ( 1  )/6'YU)  ♦.b*U«6T*  ( 1 -COS  (3. 1 4  1 

_ I _ Ylll/OtL(J)  )  1 _ 

CALL  UL 

_ UA(l>«UAC(l)>Af(l*.2«Ut.C(l)»»2/AL»02>/At: _ 

UA(2)«UAC(2) *Al*(l*.2*0tC(2>**2/At**2) /At 

L _ PRIWT  INITIAL  l/MLntt; _ _ 

WRITE  (6t?0U2)  ULCd  >  tUtC  ( 1 )  fOXC  (1)  tOFL  ( 1 )  f  UE  ( 1 )  tOA  ( 1 ) 

_ T&«0. _ 

WRITE (St^uOEl 

_ WHITE  (fei20U3)  Ui  iPI  .  Al  iKnJl  «  I  1  .UtC  ( 1  >  .HL  .  At  .Knot  .  TE  .  TK.Pt;  .  Aa«WHO  1  1 

1  ) »TTflK.PS»*S.KHOS»TS 

_ wRlTE(b.?0U4)  utm.Tn.rli  I* _ 

WRITE (6»?OOA1 

_ WHITE  (ht2(jQ71  (  (T  11)  tU(l)  t  YT  (II  lUCIll  .RHOm  >  .lal  ,m1 _ 

WRITE (btPOOS) 

_ CALL  RRlNT(-iiL3) _ ^ _ 

C  CALL  EaECOTIVE  SUHKOUTlmt 

- 3  CALL  loop _ _ 

2001  FURRaT (//lb* '  ITtRATIO^S  UbtTs • *F 1 0 . 4 t •  UT AU^ • «F 1 0 .6/ 1 

2002  FORMATl/lhAt' - UElIA _ Ut _ uEA»/.»  cqmhhessihle  «.Fln.b. 

1  2F10.2t/'  iHCOMPRESSlbLE*  «M0.6t2F  10.2) 

2QQ3  FOHRftTt* _ U _ e. _ fi _ dmi _ Tt 

1  /  •  INF  •  .3F  10. 2.F  10. 7tH0. 2 

- 2 - C-1 _ ED6Et.3H0.2.F10.7.FlQ.2 _ _ 

3  /  •  MALL'i3F10.2tFin.7tF10.2 

- 4 - 2_! - &TA6».3H0.2.F10.7«hlQ.2  1 _ 

2004  fORMATI/*  OEL*  TmETA  N  CFV»4fl0.b> 

2005  PURMATIM - K _ »_C _ U_L _ OU/OA  C  OEL  C  PEL*  TH  h _ 

ICF  NU  X  Ut  1  UU/OA  I  OEL  1  UTAU  Ub^T  TW'/) 

-^006  fUHMM  1* - 1 _ LI _ U _ lie _ RHu  t  > _ 

2007  format  12 (F 10.6 *F 10.2) .7  10.7) 

20QH  FQRMAT<////1 _ 

STOP 

Enu 
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_ SUMROUTINt:  INHUT - - - 

C  INPUT  KEAUS  ALL  IhPUT  tACtPI  IHt  UTAU  OK  UELTA.  TmE  VELOCITY 

C  nFHlWATTVES  ANU  X  toRIU  ARE  SF  T  UK  Ai^U  PWlNT  InTEHVAL  OECIOEU. - 

COMMON  YTn01)*Y(101)*UC(101>fU(101)  tRnO  <10nfXT(201)fX(?01)» 

_ 1  uEC  1201  >  .UE  (20n  .UXC(201  )  «UX(2Qn  .OLC  (2Q11  lUELLZOl  1  tUT  AU 

2  UTPJZOl)  »K<*tX6»xr  tUIfPl  »  Ai  tKtiUl  »PE*  At*KHOE»IM»M»C»OX*OY»TH» 

_ n  nPI  n.Ta.llHFT.KlMLi.E  <A>  . tJM. Kri  .  f  T  .  T  I  .  Ttn H _ : - 

COMMON  /mmm/  LS«L0*LT 

_ HEAD  l&.aon  M.N _ _ _ _ _ 

HEA0(5tb01)LS«L1 tLU 

_ HEAD  <S.»U2»  C.Kl.HriOl.^Jl.ULCtl)  .UM.XUtUTAUIJJ - ^ - 

REAO(Stb03) (UEC(i) tI«liM> 

C  SETS  PHiKf  INTEmWAl— K4  KOH  X  AND  Kfc  EQH  Y - 

K4al 

Kti»N/15 - - - 

0EL(1)>0LC<1) 

_ DVULCt  1  _ _ _ _ _ 

DX>C/<M-1) 

_ XTtll»AO _ _ _ : - 

XT(M>»XO*C 

XCllbXQ - - - 

UXC(l)»<-1.5«UEC(n^2.*UEC(2)-.b*Jt(;(3)  )/0X 

Mi«M-l - — - 

00  1  J>2tMl 

_ XT(J>»XT<J-1)*0X - - - 

I  UXC(J>»(UEC(d*l)-UEC(J-l) )/<2.*0X) 

_ llXCIMl«(1.«i*UEC<Wl-2.«UEC(M-n».&»UEClH-2»)/QX  - 

UEC ( M* 1 ) *2 . *UEC <M ) -UtC ( M- n 

aol  FtJRMAT(31&) - - - ^ - 

002  format (0F1O.U> 

a03  FORMAT  t/F  10.01 - - - 

RETURN 

EnO 
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_ SUbHOUIlNh  LOOP _ 

C  EXECUTIVE  SUBPOuTInE.  CALLS  NAUV(AUWANCt  KlTH  VELOCITY  SPECIFIED) 

-C - iiAOV  lAPVAtiCt  UlTH  ^TALI  SPECIFItJ)  ANU  OAOV  (advance  iilTH  OEI  TA _ 

C  SPECIFIED).  ALSO  kEADS  UTAiI  OK  DELTA  AND  CALLS  INTERPOLATION 

-£ - SlIPHflUIlNtt _ EXIT  IS  EHQM  MAIN  ON  HETUHN. _ 

COMMON  YT(10l)»Y(10n»UC(101)iU(101)  .MHO  (101)  .XT  (20 1 )  •  X  (20 1  i  « 

- 1  UeC(201)tUE(2Ql)tUXC(2Ul).UX(2Ul).DLC(2Ol).DE(  (;»01)  .urAU(?OH  . 

2  OTP (201) *X.K1 »AA*K6.MtOI tPlf A1 tPhOl.PE.AE.PHOE.N.M.C.OX.OY.TH. 

- i  DELD.T«.Ut>El.bNU.E(»>.UM.HH.rT.Tl.Tt.M _ 

COMMON  /MMM/  LS.lD.lt 

_ 1  CALL  NADVIL?) _ _ _ 

2 

_ _ _ 

X(K)»X(K-I)*OX*Pt*AE/Pi/Al 

- If  (X.(»E.M  ICALL  KKINT(1.&3)  _ _ 

IF  (Kl.6E.K4)  Call  PKl  nT  ( - 1  .  >.4 ) 

_ 4  If  (K.EU.LS)  60  10  S _ _ _ _ _ 

Call  nadviep) 

- S  If  (LS.EU.LTtbO  TO  bfi _ 

C  UTAU  INPDT 

- BEAD(5.1QI1  (UlAUtJI  .j.LS.LT.yi _ _ _ 

CALL  LA6R4(UTAU«LT.LSI 

_ LISLT-I _ _ _ _ _ _ 

LN«LSaI 

_ PO  1»  I»LN.LT _ _ _ 

16  UTAD(l>«(UTAD(l*l>«2.*UTao(l)«ijlAiJ(l-n)/4. 

_ LI»LT-i _ _ _ _ _ 

CALL  SADV 

_ XN1TE(6.102) _ _ 

GO  TO  1 

£ _ delta  input _ _ _ 

50  READ (S. 101 > (DEL (ol . J«LS.L0.2) 

- CALL  LAGR4(DEL.LD.LS> _ _ _ _ 

LOaLU'l 

- DO  51  I«LS.I  U _ _ _ 

51  0EL(I)«(DEL(lAl)42.*DEL(I>ADtL(J-)) )/4. 

_ LU«LD-1 _ 

Call  oadv 

_ WmTE(fe.lU.>» _ _ _ 

60  To  I 

_ 101  FOMMAT(/F10.0) _ _ _ 

102  fohmatc  velocity  specified*) 

_ 3  return _ _ _ 

END 
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_ SUMROUTtNl:  ufcxi'i  lDtLm.UlftH.Ut.K’tK9.1I) _ 

C  FOt^MS  ICItMS  FOH  Ubt  IN  THt  WUNlit-KUTTA  SUBROUTINtS. 

_ COMMON  YTiini).vnoi>.iiC(ioii.unQi).MMUHQn.»T<?Qii.x(?ni>. 

1  UEC  (2U1)  tUE  I2U1 )  «UAC(201)  «UX(2U1)  tOLC(201  >  *UE.L<201)  fOTAU<201)  t 

_ ?  UTP(2un  .H.M.X«>.Xfe.N/.Ui.Hl.Al.aMOf.PE.AF.RMOfc.«N.MiCiQX.DY.TH. 

3  UfcL0*T«.UHFT»bNU»E(4) .UM.HH* T T f TI * Tt»M* TS 
_ COMMON  /DiJ/  ULPMlPnil _ 

yb»o, 

_ d&«0.  —  . . 

B7a0. 

_ f  JlElLi _ ; _ 

Fb*0. 

_ ^  tsQ. _ 

AllaO. 

_ Alg«0.  _ 

A13S0. 

_ A21c[j. _ 

42230. 

_ A233U. _ 

S3TT/TS 

_ L>«DE.La/(N-n _ 

ULQ>UELu«Ab$ (U T AM) /Gnu 

_ nu3)-l(i;m) _ 

H20SO. 

_ F203(2.S/(1»0LU) )-d2Q _ 

Ft>03F104UL(.i*FcO 

_ A13L3UEH«Utm«S _ 

A23L  =  A13L*l<tCw/4r 

_ H143 _ -UTAk«ab&  (in  AM) _ 

U243bl4 

_ PlE»3.l4l!»«i> _ 

UbETsUEM-uTAM*F )u 

_ UTQ3UTAm«ABS  (UTam)  /Gn'J _ 

DO  1  I32»N 

_ Y(1)«T(1-1)4U _ _ 

Yi.U*r  (  1  )  *A[>S<UTAr')  /Gnu 

_ U(I)3UTAF«F  1  (rLb.)4.S*Ur.ET*<l-CUi»lPlb>Y(l)/QbLU)) _ 

AMU»J*(-1)**1 

_ IMl.EO.N)  AMU31. _ 

BGA>F6 

_ B6A»F7 _ _ _ _ 

B7A3F3Y 

_ lUJ«llll«ii.aSigiAtLL/&JU_ _ 

F1Y3H  (YLu) 

_ fijgOEQ _ 

1F(YLN,1  T.bO)  e«:0s-(1.2b4*YLO-l.S02)*tXP(-.37*YLQ) 

_ F2Y»(2.B/tl*YL0)  )-r«2fi _ ; _ 

F3Y*.5*ll-C0S(Plt*V(I)/UtL0) ) 

f4Y  3  tU-JLil  N .( Hlfe  *JLLI  )/ufc.LO>/2. _ 

FSY3|-]Y4YLw*F2Y 

_ Fb3F!aY-(-  3Y*FbU _ 

^  /3-06E  T«Y  ( I )  *F  A  r/OFLN**2-UT(4*F  3  Y*F  20 

_ UY»UTf)»F?Y4Ubt.t«F4Y/0EL0 _ 

«b»bb4  (b‘jA*Fb)*lJ/? 

_ bb«Hb4(bbA*F  n«o/2 _ 

b73M  /4 (b7A4F  3Y) *u/? 

_ All3All»(U{I)>Fb-UY»t<b)»A''iU _ 

A21*A2l4  (U  (  1 )  *Fb-UY*B‘j)  *A’iU  *Y  i  I ) 

_ A123A124 (U ( I )*F  y-UY«bb)»AwU _ 

A22»a22*)u(1)*F  (-1(Y*Mb)«amU  *7(1) 

_ A133A13* (U( 1)»F3I-UY»b7)«amu _ 

Ad33A234 (u ( I ) *F  3 Y-UY*b  7 ) *AMU  *Y ( 1 ) 

H£T3  HtTA(1.0tLN»UTAK.uEH»II) 
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1  *ij  r 

J  _  _ _ 

=  _ _ 

_ _ _ 

h^4  =  tl/?4*l>/  3 

A i 3=a1 3  *DZJ - A ! 3  L  _ 

A^isA^^jAO/^-Ac  .*l 

AJfi?AZl/'Jtl..',L_  -.  _ _ 

Ai:c»A»»r/('h  t  ij 


_._aiii=AZi/12tLU - 

l  O 

_ Ih_  Ll.JlL_  MV  Jni-d - 

It  .<)  toO  10  3 

_  OAl>»U*  IfN*  1  )  -UA  (A  I _ 

!»■  (ii.to.n  jxu=n 

LlXO^bX  liv)  ♦UXO/t  lii  ) _ 

rtl4s-Al3*|.iXI)*hl<f 

=  A^i^iJAU _ 

UtM*l.iXu 

_ ntLQg  (  Ajjl  1  >  /  (Al^»Af;l-Ali«Ai:!^) 


jliJJiiiN. 


UTl)»UTH  («  ♦  I  )  -UTk  (►  I 

J±_ XliitOi  J 1  uTui  =  n 


IJTHPaUTH  IM  ♦oro/t  (in 
Ml4«ril4-A 1 1«U IPf _ 


M<;4*rii?4-A^l*UlPP 
UTAP«OTkP 


OtL(J=  («14*A>»J-d«!4*rtl  .3)  /  (A  l2*4i?J-Ar'<:*Al3) 
Ot.P»(Hl4«Ai;?-h^4  »Ai;*l/(A13gAf’?-A23«Ali>l 
><tTUKN 


II-dl.LU.l)  0LU3l>. 
ULHP8!)LHW(K)  ♦OLU/t  111) 
Hl•^»ISl4-Al?*UL^'^' 
bi>4a0ij4-Ai^^«0LPP _ 


OtLQsliLHP 

urAP«  ia^3«hi4-aij«p?4)  /  (aiiaa;?j-ai3»a2H 
UfcPs ( A21*p14-a1 1»K24) / (A1 1*A23) 

KhTUt^  * _ _ 

tNU 
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_ SUP><UUTlW|r  NrtUv(<) _ _ _ 

C  KlIHOA  MiCM  AOv*>N(,c.  tXIlrC'dXL  VtLUCilr  IS  S«>tCUItU.  CALLS 

£ _ l.'fcS.lXJ'.'4t<._LykKtI L H:.v lb  I^b  utHivAiivLS. _ 

COMMON  TT  (  lun  »T  U  01  )  tOL  (  1  P  J  )  »U«  lou  »KhO(101  )  t  A  T  ( ifU  ]  )  *  A  ( <rU  1  >  « 

- 1.  -U£C-L^i)jJjfc  .oAClgPi)  tOAlfiOi)  tOLC(ir>1U  «Dt.L(i!01)  tJTAUI^Oll. 

li  UrP(?Ul)«A>AI«AA.A^.K/,Ul*MX,Al«>4riOl  t  M  t  »  At  *  MHOk  •  N  •  M  t  C  *  OX  •  OV  *  I  H  t 
_ 3  F  »oNOtL  <  A)  tOM»KWt  l  T  t  fl » l£tM _ 

_  _ 

UTAUMsu. 

_ J)Ul  =  0j _ . _ _ 

_ U»OXAPt:»AL/  _ 

C  RONbt-AOTIA  AOVanLL 

_ KO  1  _ 

UtLPsOtL  (A  )  ♦ULLI-*(j/t  ll ) 

_ UlAUPaUl  AU(K)*UFmIIP«»U/|-  (11 _ _ 

UE.H*UE  (A)  ♦UF  (  1 ) 

_ CiOJ.  -ULhlli  (QtJ. P.oTAUPtUtLP.liIl _ 

OOlsLII>i*A.  (  I  )  ALtLM 

_ 1  iJUl«UUi*ACl)*OIALiM _ 

J»A*1 

_ DLL ( J) =UtL ( K ) ♦OOIAH/A _ 

UTAU(J)=U1AU(A)  ♦Ulll*0/o 

_ CALL  EaT<J> _ 

OC.(J)sUbC(J)*»l/<«F 

_ UA(J)«UXC(J)«AlA<l*.2«U>-C<J)«»;^/ALAA^l/At. _ 

!► < j.Gt.Mifao  ro  j 

_ L5J±J _ _ _ 

3  IU<L)»UaC(L)*aI*  (l^.^AutCdlAAii/AtAAifl/Afc 

_ PtlOHN  i _ 

tNO 
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_ SUriHOuHMt  DAuv _ _ _ _ 

C  HLINGA-MJTH  flliVANCt  H(HJNU«HY  LAYtK  TMlCKiXtSS  IS  SKtCIFltU. 

- CUflwtfN  YT  tion  It  tlQlI  tijC(lDi)  lUdOl)  iKhOdOl)  .XT(?0  1 )  .X  typi  I . 

I  uec  (?on  *ut  (iJO)  )  tO;iC(/?Ol )  ,UA  (2U1 )  tiJLC(20I  >  tUtL  (201)  »UTAU(201 )  . 

- a  UTP  t  A1  tHnUl  tPt.  ALikHQfc.tN.M.C.OX.I)Y.Trt. 

3  UE.LU«  r  Hi  tUatT  .GNUtt  (<*)  nJMfWWtrTtTItTEfM 

- COMMON  LStLUtLT _ ^ _ 

COMMON  /l)U/  i>LMKi201) 

- «<HlTiL(b.«<Ul) _ _ _ _ 

J»K 

- UtfLTgUfc(J)-UTAU(Jl»Fl  lAbSiUTAUIJl  )  <H>tLL  I  J)/GN1]>  _ 

CALL  PH1NT(-I*fc«) 

- B  0LPK(K)»(^t■L(^♦l)-Ot:L(^l  )/(L>A»Pfc*AL/(MI*An  1 _ 

OLPH  (K*n  a  (UtL  (n*2)  -OtL  (K*!  )  )  /  (l)X*ML*AL/  (Ml*Al )  ) 

- OQ  2  KK.LS.I  0 _ _ _ 

K*KK 

_ 0«0X*PL»AL/(PI«An _ 

ULLQso . 

- mi:p»u. _ , _ _ 

UTAOP=0. 

- GTlai.. _ _ _ 

UUlaO. 

_ 00  b  I  a  1.4 _ _ 

OLLOsOtL  (n  )  ♦UtLO*i)/fc  (  1 ) 

- : - UrAUPaLrAmK>»UTAIlP»0/ti  1  1  1 _ 

UtPaoe (A) ♦OEM*U/t ( I ) 

- CALL  OLMlV  (QLm.uTAUPiuLHi  J.l) _ 

UIIauTI*t<I)*UTAUP 

- S  UUlaUOltL'd  )*OLh  _ 

J»R*1 

- UAlJ)aUol/h. _ _ _ _ 

lJTA0(J)al)rAU(K>*uTI*0/b. 

- UL(JlaOL(M  ♦OUI»u/b.  _ 

CALL  EXTIJ) 

- ULPH ( Jla  fuLl  (Jtl )-ULL  IJ) ) / < OX >Ph » AF / 1 P t « A  I  1  1 _ 

IE (J.Gt.M)UO  10  J 

- _ _ _ _ 

3  OLPP(L)»(OtL  <L*l)-OhL(L) ) / ( Ox*Pt* At/ ( P I • A 1 ) ) 

- Lit-HJlauK  (J)«At/i.r _ _ _ 

Kaj 

_ KlaKl*l _ _ 

X (K) ax (X-1 ) *OX*Pt«At/Pi/A I 

- IK  (M.bE.KAJ  CAM  PKli^H-l  .A?) _ _ _ _ 

t  CONllNUt 

_ HtTOHN _ 

‘*01  FOWMaT(//*  StP*.rtATtO  MtOlON  utLTA  bPtCIFltO*/) 

- LJXU _ 
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_ SUtfHOUTlMt  _ 

C  HUNGt-KUTlA  AOVuNtt  KHtN  fKlCFlON  VtLOCIIY  IS  SPtCIHfcO. 

- CUMMGN  YTHOH  tY  tlQl)  tULdOn  .  U  1 1  U  1  )  .  WHO  I  1  0  H  .  XT  1?0  H  .  X  ( ?U  1>  . 

1  Ut.C(?01  )  (Ut  (<!0U  «UXC(^01)  .UX  (2U1  I  tOLC(?Ui)  tl>rL(?01)  *UT*U(201)  • 

- 2_UTPli;QJJ  iN»K_i..RH.KbtKdUl.f»itAltHftQl.Pfc.AE.HHOt.«N.M.C.QX.tJYtTH. 

3  l)liLD«  1 1^  t  UmE  T  »  Gr*0  *  1 14)  «UMtHr(tlTtTl«Tfe.<H 

_ COMMON  /KMM/  LS.LOilT _ 

\XHI1  rt. 

_ JSi _ 

nHtT=iJt(j)-urAU(j)*e  i  (/ibS(utai/(J)  )*otL<J)/GNU) 

_ CALL  -gia  _ 

«  UTM(K)  *(IJTttU^^  ♦D-UTmUIK  )  )/(l»X»Pfc*Afc/(Hl*AI)  ) 

_ UTP(^♦l)»(UTAU(^♦?l-nrAU(K♦ll  l/tOX«>*t.«Afc/tHI»All  1 _ 

1)0  2  KK«LS*L1 

_ IL»»M _ 

UaOX*Pt*At/ <P1*A1 ) 

_ Ot-LPaQ. _ : _ 

UfcP«0. 

_ UTAUPab. _ 

001 aa, 

_ UUlBiU _ 

OO  S  1  a 1  •  4 


UTAIJPaOl  AU(K  J  (  I  ) 

_ UtPaU^L^^)  ♦UtlPOO/tm _ 

Call  otNi V (Ofc.LP*uTAOPfOtPt-l • l ) 

_ DOlaOUl^fc  li  )  *L>bLri _ 

b  UUlaUUi*):  ( 1 )  *UtK 

_ _ 

UA<J»  aOuI/t. 

_ OliH  Jl  aof.LlHl  ♦ULIAU/O. _ 

Ut(J)  *Ut  (N)  ♦UUl*U/‘»» 

_ Call  KXI(J) _ 

ltTP(J>a(UTmj(J*l  J-Ut  AOIJ)  )/(l)X*PL*At./(»-I«Al)  ) 

_ Ufe.C(J)au(-  1 J  )  »At/Ar _ 

le (J.GC.M)GO  10  3 

. . L*J±1 _ _ ^ _ _ 

3  urp(L>«(UlAli(L*l  J-')TAU(L)  )  /  (DX^PfcAAt/ (P1»a1  »  » 

_ K«J _ 

KlaKl*! 

IF  (Kl.UL.KA)  CALL  PP1nT(-1«N2) 

_ i  CmUJYUJL  _ _ 

Pfc I Ukn 

_ 9fll  FONMAft//  ■  StLPAMATLU  PLGIOn  UTAU  SPLCIULm/l 

ENO 
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SUD>«OUTiME  BL 

.£ - CQ^lf^UTtS  Ul&fcLACfc^tM  AiSl>  ThlCANitLSSiii.. SHAPE  ^A^.TOB  ANU 

C  SKIN  I^BICTION 

- CUHr^UN  YT  1101)  t  YIIOI I  .uciioi)  .111  nn  1  .urtu  ( i  o  1 1  .xT  (>»n  n  .  a  tpn  i  > . 

I  UtC(ifOl)  fUf-  (<IOn  .'JXC  «i!0l »  »IJA  (201)  *LiLC(2ni)  .DtL<?01 )  »UTAH<20n  , 

- 2  UT>>(2(>iI.iRt'SltK‘tiKtH<SyjUl-tKliAiiHrimiPttAE«kHOH.N.M.C.UX.tjY.TH. 

■i  DELL)*  iMtUBE  T  *(;r*u*k  ( 4  )  *Um*HM  *  T  T  «  T  I  *  Tt  tH 
- L»L1.U«DY _ _ _ 

TrisO. 

- DP  7  _ _ _ _ 

23«1  .  -  (WHO  ( I )  *00  ( I )  )  /  (wHO(:*UtC  (K)  ) 

- 21»Z3*(WWU(l)*UL(n  )/(r^MOL>JfcC(^i  1 _ _ _ 

Z3«?3*(3*  (-l)**n 

If-  tl.bWiN) _ ^3.Z3/4 _ _ _ 

(3* (-1 )**i ) 

- li— tlitWiNI _ /l»Zl/4 _ _ _ 

DELD*()tLn*73*uY 

_ 7  TM»Tfi*7l»0Y _ _ _ 

DtLUsUfeLO/3 

— — : — iti«Tn/3 _ _ _ 

M»UELD/1M 

- lM-«UTAO(M*Ad:»(UlAIJ(K)  _ _ _ 

WtTUwN 

END 
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_ _ t-UNCriuiN  MfcTA(j<up j ) 

C'  Co’i^HUTt-.S  tljUY  ViStUSiTY 


_ _ CUtf'li^jY.Jr.UlJ)  li  tY(ioi)  .ucdOl)  ,>i(  101)  «Kno(ioi)  ,aT(^01)  .x(^Qn  « 

1  UEC(?01)  >ut  <^01 )  tOXCC^U  1.)  «OA  (/^Ul  )  tULC  (^ai)  fOLL(^Ol)  «tJTAU<?0n  « 

3  UELO«t'iitihEli  »t  (ft)  tiM»WK*T  I  »TItTt»rt 

_ IE...  (UH«lL«U) _ _ 

Ir  (J.bi.^*)  «o  10 

_ U«UP/ (b-1) _ 

TMaO, 

_ _ _ 

l»U  3  lSii»N 

_ TM8TM«(l.-ii(l)/m,)»(3*(-l  i««l)»i.j(i)/mi) _ 

J  l»ELl)=l)tUl'*  ( )  .-IM  1 )  /i.'U)  •  t  3*  (  -  U  **  1 ) 

_ Tm«Th«U/3. _ 

t)tLU=lJtl.L)*Li/3 

_ »Jj-3 _ 

TiiJ*Ok-*AMS  (UP )  ♦Hhyr 

_ eA?L4) _ 

IE (Ua(a) .Lf.U.) 

_ lPA»-KMUll»JU»tJX  IM _ 

f'l»ArtS(OELU*H*/  (  IW^IS.)  ) 

_ -It  (  tLL.LWiAl _ 

lPKs.Ul3*.U03H»tAK(-(L)LLU/T»()*>'X/lb) 

_ !LL«--J _ 

2  lE(X/.bT.O)  bu  TO  <* 

_ Tl«U(J)/AhS(UK) _ 

IE  (Tl.bl  .(?b.  )  bu  TO  4 

_ T2»I-*.0i>33«<tAE{.41»ri)-<l».41411*.b«(.414T114»2)  ) _ 

T3«Pk*  ( 1  ♦b.b*  I  Y  <u»  /UP)  •*(!.)•*  <-l )  *l)u4utLU/bNlJ*  1 

_ IE  (T3.Lt-.T/)  A7ai _ 

li^fcTAaA^*lM  (T2tl3) 

_ wtlurti.' _ 

4  HtTA=PX»(l*b.b*(Y  (J)/uP)**b>**(-l)*UU*Ut,LO/(iNU  ♦! 

_ HEfUU.Y _ 

1  IE  (j.bl./)  bo  in  h 

_ f)0  b  I»2«N _ 

Hal 

_ IE  (Udll  .bE.Q.)  uQ  lU  f _ 

5  CONTlNUt 

_ 7  OELOaO _ ' 

OaUP/(N-l) 

_ UQ  b  lall.N _ 

8  ueLl>*0tL04(  d-(L>d)/00)  )*d-(Ud-l)/UU)  )  ) 

_ UELUaOfc.Lr)a.Ol  jauu»U».i»/b.sU _ 

h  RET A3ULL0*  d  ♦:».*>•(  Y  (J) /UP)  4*6 )*♦(-!) 

_ return _ 

Enu 
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_ SUBROUTINE  PHINT  IKg.«l _ 

C  COMPUTES  compressible  VELOCITY  PHOFiLtS  ANU  PRINTS  OUTPUT 

_ common  YT(lQl)  tYUOl)  lUCdUl)  .UtlOl)  .RMO(lQl)  .XT(2Q1)  .A(gOl)  . 

1  UEC(2UI ) «UE (201 ) lUXC (201) tUX (201) • OLC ( 20 1 ) « DEL ( 20 1 ) fUTAU(201) « 

- 2  UIPt2Dil  t^^«Ki,^4nK6.K7^UItPI^AItRMaltPE«AE.RHQE.N.M.C.OX.nY.^ri. 

3  0ELU«Tw«UbET<(>NU«E(4)  *UM,PNtTT«TI«TEtM 

_ dimension _ Ldol.pol) _ ^ _ • 

DIMENSION  K0U(2<») 

_ R0tZC,ZEl.l4.1/B«</b/AE)««2«Kl-(7C/7El»#2> _ 

IF  (K2)  ltl«3 

_ 3  M»1QQ _ _ 

1  DY«DtL(R)/(N-l) 

_ U^ET«UE(Kl-lJTAU(^)^Fl  (  AhS  ( LJT  At)  ( M  *UEL  (K ) /SNUI  ) _ 

RM0<1)>RH0E/R0(U«UEC(K) ) 

_ DO  *  1*2. N _ _ _ ^ _ 

RMO(l)»HMO£/HO(UL(l) »UtC(K) ) 

_ 4  YT(11«YT (l-l)*Ai«RH0I»UY>(((l/Rrl0111 >*(1/RMQ(1-11 ) 1/AE».S) _ 

OLC(K)»YT (N) 

_ DY«0LC(H)/1N-1) _ _ 

DO  6  I«2«N 

_ YIll«Y(i-lJjt  IHMliU)  tHM0(l-l)l»AE«UY/(RM0I*Al«21 _ 

U(I)  aUTAU(K)*P  1  (aBS(UTAU(R)  )  *Y  ( I ) /(jNU)  ♦  . b»UB£T*  ( 1 -COS  ( 3. 1 4 1 59« 

_ 1 _ Y(I1/DEL(K)  )1 _ : _ 

6  UC(I)BU(I)*At/AI 

_ CALL  BL _ _ 

00  8  lBl»NtK6 

- a  UliMM  ULUy  UtlKl _ *1000 _ 

L(N*K)«V99 

- BHlTt(6t2UQl)MiAT(M  lUEClM  tUXC(M  tULC(M .DELD. 1 M.ri. Tw.GNU. X (K > . 

1  UE (K) tUX IK) t DLL (K) tUTAUIK) •UUtTtTT 

_ IF  IKl.bT.50)  SO  TO  2 _ 

Kl»0 

RE-TUKN  1 _ 

2  CONTINUE 

_ IlHil _ : _ 

UO  14  Isl.N.Kb 

_ 

14  K0U(I1>>I 

_ MrtlTLtb.2002)  IKQull)  .l«I.N.Kfal.(lI.(L(l«II).I«l.N.K61.Il»l.M.KA> 

2001  format  (i4,F7.4tFb.0tF8.O.3Fa.S«r-S.2«2F9.6tF6.3*F6.0»Fe.0* 

- 1 _ Fa.5tF7.3«F/.ItFS.O) _ 

2002  FORmaT(»1  velocity  PH0FlLES»//»5Xt21I5/210 

_ 1 _ L22I&/11 _ 

RETURN  1 

_ EJiU _ _ 
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_ SUBHbUTlNt  L*bH4tAtHHL) _ 

C  FOJHTH  OkUFk  LAbR/tivtilAN  li^TtRHOL AT X ON 

_ UlMFNSiUN  XXXiOIl _ 

L»L*3 

_ _ _ 

00  1  J«L«N*? 

L«L-J» 

X(L)»(AtL*l)»A(L-n)/2. - 

X(M-l)»<*<M)»X<M-2) )/2, 

_ L»L-i _ 

RtTURN 

£nQ _ 


_ SUBROUTlNt.  EXUJi - 

C  COMMUTES  toot  CONOXTIONS  ANO  iNCONRHtSSIBLE  VELOCITY  PKOFILES 
_ COMMON  YT(IOX)  .Y  HOI)  .UC(XOH  .0(101)  iRHOdOll  iXT  (?011  «XI20iy  1 _ 

1  UEC(20X) «Ut (201) «OXC(201) tUX(20X>  tOLC(201) tOEL(?On lUTAUlPOl) « 

2  0TP(2U1)  ♦K«M«^^.K6«^T^Ul^^»^^AXtMhQl^HE^At.^ilM0t»^■MtC■0X.DYlTri. 

3  OELO«  T  w  «  UBt  T  <  GNu  « t (4)tUM*HR*TT«TI«TE«M 

_ MQ(/C.2t<«l«W7B«l2t/Atl«*g*H-(2CyZtlM2X - ! - 

6AMal.4 

_ GAM1»6AM-X _ ^ - 

AtaSCiRT  (.?*  (UX**c-UEC(J)**?)  ♦AX**2) 

_ Rt1QEa(HM01«»OrtMl*6AMl«(KH0ia*SAM/PIia(OX**?»UECtJl»*21/tg*fiAMna» 

1  (1/uAMl) 

_ PtaRMOfAf  «2/1.H - 

Tt*RR*Pt/PHnt 

_ ONUaUH»TE/T I/MMUt - 

RHU(l)aHHUE/HO(0.<OEC(JI ) 

_ TTaPE»HH/RMO(  H _ 

RETURN 

_ &iiU _ 


_ t UNCTION  EHZl _ 

C  LOGRlTHMlL  function  FOR  VtLOClIY 

_ Li _ a;,;.S«AL0U(l*71*ia.l _ 

IF (Z.6T.A0.)  00  10  1 

_ Harl  -fJ.3»<»7»^.l>»EXPt-. 37*21 

1  RETURN 
END 
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Aij 

Cf 

q. 

Fi 

H 

P 

S 

T 

U,V 

u.v 

U0 


Ut 


u+ 
X.Y 
X,  y 
y* 

P 

5 

«• 


E 

6 


NOMENCLATURE 

Sound  speed 

Coefficients  (Eqs.  19  and  20) 
Skin-friction  coefficient 
Pressure  coefficient 
Defined  functions  (Appendix  B) 

Shape  factor 
Pressure 

Enthalpy  parameter 
Temperature 
Incompressible  velocities 
Compressible  velocities 
Wake  velocity 
Friction  velocity 
ulxij. 

Incompressible  flow  coordinates 
Compressible  flow  coordinates 
|ur|y/i; 

1  +  elv 

Boundary-layer  thickness 
Displacement  thickness 
Eddy  viscosity 
Momentum  thickness 
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V  Kinematic  viscosity 

p  Density 

Tw  Shear  stress  at  wall 

SUBSCRIPTS 

AW  Adiabatic  wall 

e  Edge  of  boundary  layer 

t  Stagnation 

<»  Free  stream 


44 


